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Symmetry restoration processes during the non-equilibrium stage of "pre- 
On ■ heating" after inflation is studied. It is shown that symmetry restoration is 

very efficient when the majority of created particles are concentrated at en- 
ergies much smaller than the temperature T in equilibrium. The strength of 
symmetry restoration measured in terms of the equivalent temperature can 
exceed T by many orders of magnitude. In some models the effect can be 
equivalent to that if the temperature of instant reheating would be close to 
the Planck scale. This can have an important impact on GUT and axion 
models. 
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In accordance to modern cosmology the outcome of the earhest stages of the Universe 
evolution, which predefines its modern appearance, is determined by fine details of the 
dynamics of a particular scalar field, which is called inflaton. For a review of inflationary 
models see Ref. [0]. The chaotic inflation is a typical model which has essential common 
features. In this model the inflationary stage itself persists till slowly decreasing inflaton 
field, is larger than the Plank scale, Mpi ^ 10^^ GeV. During this stage the Universe 
expands exponentially and the room for the future matter is created. This stage ends when 
the inflaton field reaches ~ Mpi and then the field starts to oscillate coherently. Coherently 
oscillating field can be considered as a collection of unstable inflaton quanta at rest, so it 
decays to all particles it has coupling with, and the matter is created. 

With the assumption of "instant" reheating the products of inflaton decay would ther- 
malize on time scale negligible compared to the rate of the Universe expansion, so the 
temperature after reheating would be T^^ = p^/f?*, where numerical factor g^^ includes total 
number of degrees of freedom and is large, (7* ^ 10^. Here is the initial energy density 
stored in inflaton oscillations ~ A0^/4 ~ AMp^. It is relatively low since the inflaton 
self-coupling constant has to be very small A ~ 10~^^ for the induced density perturbations 
to satisfy observational constraints. In reality, the reheating temperature will be even much 
smaller since the reheating is not instant and while the particles thermalize, the Universe 
expands and cools. 

The magnitude of the reheating temperature after inflation is considered to be important 
as, for example, this will determine whether or not certain scenarios of baryogenesis in Grand 
Unified Theories (GUT) of strong and elecktroweak interactions will be successful which 
requires the reheating to be up to the GUT scale Mx ~ 10^^ GeV. Another important 
issue is the occurrence of phase transitions. If the reheating temperature is larger than 
the Grand Unification symmetry breaking scale, then the corresponding symmetry will be 
restored (for reviews of phase transitions in GUT, see, e.g. "^^^ subsequent cooling 

will be accompanied by a symmetry breaking phase transition which will proceed in different 
horizon volumes independently, resulting in creation of topological defects: domain walls, 
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strings and magnetic monopoles. This means that the problem of monopoles and domain 
walls 1^ will be resurrected, ruling out the corresponding models. 

The question of symmetry restoration is interesting not only in connection to topological 
defects. Another important aspect is whether the Peccei-Quinn (PQ) symmetry is re- 
stored after inflation, or not. PQ symmetry was introduced to explain the apparent smallness 
of CP-violation in QCD. Pseudo-Nambu-Goldstone boson resulting from the spontaneous 
breaking of this symmetry, the invisible axion, is among the best motivated candidates for 
cosmic dark matter. The combination of cosmological and astrophysical considerations re- 
strict the relevant symmetry breaking scale, or axion decay constant fa, to be in the narrow 
window 10^° GeV < fa < lO^^GeV (for a review see |^, note that the upper bound on /„ 
does not apply in certain inflationary models If the PQ symmetry is restored after 

inflation, then the axion field will not be constant throughout the Universe, but will have 
independent values in different horizons. These fluctuations in the axion field are trans- 
formed into density fluctuations of order unity at the crucial epoch when the axion mass 
switches on at T 1 GeV, leading to the existence of very dense axion miniclusters P,pi] 



which may be observable [|rT[]. This also shifts the main source of axions from a coherent 
misalignment angle to decaying axion strings [Q. 

Traditionally, the answer to all of the above questions was associated with the value of 
the reheating temperature after inflation. The purpose of the present Letter is to show that 
the reheating temperature is actually irrelevant here and all processes of interest are even 
more efficient while the system is still out of equilibrium. This is more or less apparent for 
the baryogenesis since one of the necessary conditions for baryon number generation is the 
absence of equilibrium. We shall show that symmetry restoration is more efficient too in a 
non-equilibrium state generated at the final stage of inflation. 

Recently it was realized that in certain cases the decay of the inflaton fleld can be a very 
fast process p!3HT5| (see also [l^), owing to the possibility of stimulated decays. This is also 



called the parametric resonance, for the general theory of it see, e.g. []T^. Sometimes the 
parametric resonance is considered as a very special type of decay of coherent fleld which 



can not be described as a decay of particles at rest. This is an incorrect conclusion - one just 



has to include stimulated processes in addition to spontaneous [jT8[ for the particles decay 
or annihilation. That is why parametric resonance exists only for Bosons in the final state. 
But even in Bose systems the parametric resonance is not necessarily always effective. For 
example, it was studied long ago for the decay of the axion field [|T^|IB[ with the negative 
conclusion that the expansion of the universe removes particles from the narrow resonance 
zone too quickly, blocking the entire process (to reach stimulated decays of axions bound 
in a gravitational well is not impossible in principle, but would require enormous densities 
of particles in this particular case ||18| , |20|1 ). This means, in particular, that if the infiaton 
potential and infiaton interactions are constructed in analogy with the axion potential, which 



is the case in the model of "natural" infiation I^T], then the parametric resonance can be 



ineffective. However, as has been demonstrated in [|T^, a successful parametric resonance 
can occur in the case of, e.g., chaotic inflation 0]. 

In the case of successful resonance almost all energy stored in the form of coherent inflaton 
oscillations is transferred almost instantaneously to radiation, but the products of inflaton 
decay are still far from equilibrium. This intermediate stage was dubbed "preheating" in 



Ref. [|T3. To simplify subsequent discussion let us flrst take the distribution function of 



created particles to be of the form 

f{p) = A6{po - E) (1) 

For the case of two particle decay E in this equation is equal to half the inflaton mass, 
for 2 — s> 2 annihilation E = m^, for the processes 4 — 2 of self-annihilation we have 
E = 2m^, etc. The main point which is crucial for the subsequent discussion is that E is 
typically smaller by many orders of magnitude than the temperature T of instant reheating. 
Equilibration will take time, meanwhile particular phase transitions can occur in the system 
when it is still far from the equilibrium and the distribution function is still given by Eq. 
(|^). Moreover, we shall show that in this case the symmetry restoration is even much more 
efficient. 



Before we proceed, let us specify the general field and particle contents of the system. 
First, there is classical inflaton field which we already introduced as (j). We denote as 
the effective inflaton mass at the Plank scale, which includes contribution from the self- 
interaction m'^{(f)) — m(j,{Oy + 3A0^. For us it is only important that the overall value of 
is fixed to m^p ~ 10^^ GeV by the observed magnitude of density perturbations. Second, there 
are products of inflaton decay. We denote them as rj and the mass of corresponding quanta is 
rrir,. The possibility of the inflaton decay into 77-quanta assumes that there is interaction of, 
say, the form grf(f)^/2. In reality, there can be many channels for the inflaton decay and the 
flnal answer is the sum over all species. Since the true content is unknown, we shall not carry 
out this summation, but implicitly assume that it has to be done. Third, there is relevant 
order parameter, the classical held The ?7-particles couple to the order parameter, so 
that their mass depends upon it, = m^($). The typical case is = m^(0)^ + 0;$^ 

with a being a product of coupling constant and some numerical factor which depends upon 
particular direction in the space of internal symmetries. In the case of simple one-component 
scalars this assumes interaction of the form q;?7^$^/2. 

Note, that since we have to carry out the summation over all channels, there will be 
terms when 77 corresponds to quanta of the or $ flelds, and even when 0, -q and $ is one 
and the same field. In some models some particular terms can be negligible. The simplest 
possibility corresponds to the only one dominant term (p — rj — ^. While we do not exclude 
this possibility, in order to keep uniformity we shall keep separate notations for the separate 
aspects of one and the same entity. For example, we shall still denote by rj the quanta of 
the field, etc. 

In the vacuum state the symmetry is broken, which can be described at the tree level 
by the potential Vq = — /x^$^/2 + A$<l>*^/4. This relates /i to the symmetry breaking scale 
via /i^ — A$$^. At non-zero density of particles there are corrections to this potential. In 
particular, m|(0) = cfVo/d^^ — 0) became positive and the symmetry is restored. To 
calculate the modified potential wc can make use of the fact that the effective potential 
is minus the pressure. This definition is physically transparent: if two phases can coexist, 
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the phase with the lower value of the effective potential will have higher pressure and the 
bubbles of this phase will eventually occupy the whole volume. 

The pressure P for an arbitrary distribution function of //-particles can be found by using 
the formula 

--/(&^/'^)' 

where Pq = + m'^{^) and for an isotropic medium we have T*-' = 5^^ P. We omit the 
subscript r] for the distribution function since rj are the only particles we have. 

While the field $ evolves in the effective potential, the number of particles does not 
change on time scales of interest, so we shall calculate assuming that 

^-/|^^<^)- 

is constant which fixes the normalization factor A in Eq. ([l|). The procedure is very simple 
and the result for the particle-dependent part of the effective potential is Vi($) = —P = 
N [m^($) — E'^]/3E. Since for the distribution function given by Eq. (|I]) the energy density 
is simply p = NE, and since with the assumption of instant reheating the energy density 
would be unchanging, it is convenient to rewrite our result in the form 

^i(<^') = ^K(<^')-f • (4) 

At = we recover the equation of state for massless particles, P = p/3. But p is a 
^-independent constant and is insignificant for us here. What we are interested in is the 
coefficient in front of m^($), which we denote B, 

This coefficient is positive, so when added to the negative mass square of the field $ in 
vacuum, leads to the restoration of the symmetry at large density of particles, V"{0) = 
—fi'^ + 2aB. The symmetry is restored at i? > p? /2a ~ (A<j)/a)$^ (in the typical case of 
positive a; negative a is also possible in models with several scalar fields, see Ref. [Q, which 
breaks the symmetry instead). 
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If we were to calculate the effective potential with an equilibrium thermal distribution 
we would obtain i?eq = T^/24, see Ref. ||22|. Roughly, in Eq. (|^), we would have in this case 
p ~ and ~ T. We can generalize the expression for B as been given by the ratio of 
particle density to the mean energy of particles for a case when the width of the distribution 
is finite, but still it is smaller than E. 

Now we can compare the effectiveness of the symmetry restoration at preheating to the 
instant reheating (actual reheating in the expanding Universe is even less effective). In 
both cases the energy density is the same and is equal to initial inflaton energy density, 
but at preheating E <^ T. The restoration of the symmetry is much more effective in the 
non-equilibrium state Eq. (^, its strength amplified by B ~ (T/£')^i?cq. 

Let us make an order of magnitude estimate for B in some possible inflationary models. 
The inflaton field strength at the end of inflation is of order Mpi, so that the energy density 
in inflaton oscillations is given by p ~ m^Mpj. Its magnitude is fixed by the magnitude of 
primordial density perturbations. To use Eq. (|^) we only need the estimation of E. We can 
consider three cases here. 

1) If ~ we obtain B ~ Mpj , which is equivalent to a would-be reheating tem- 
perature up to the Plank scale. Note, however, that we had not calculated the numerical 
coefficient here, which can be rather small. One example when this regime can be valid 
corresponds to the self-annihilation of the inflaton field, i.e. rj = (p dominated case. De- 
cay corresponding to the 4 — 2 processes is unsuppressed in the expanding universe if the 
inflaton mass is dominated by the self-interaction, and E ^ 2m^{(f)) ||13|| . 

This effect will cause the inflaton field to "roll back" to some extent in the new inflation- 
ary scenario (the inflaton can not roll all the way back to the origin if was the initial 
condition), which was observed in the detailed numerical simulations of Ref. ||T5|. 

2) Let us consider the axion case (or to this extent any model with sufficiently low 
value of symmetry breaking scale). The Peccei-Quinn symmetry breaking scale fa ^ 10^^ 
GeV is comparable to the inflaton mass. Therefore the contribution to the mass of the rj 
particles due to interaction with PQ field is smaller than the inflaton mass and no additional 
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kinematical constraints appear. However, the inflaton will decay into lowest zones where 
E ~ only if (7 ~ A. This happens at (7 ^ A too, but then the dominant process is 
self-annihilation of the inflaton to its own quanta. For the case g ^ X the mean energy 
of created particles was found in Ref. ||13|] to be given by E'^ ~ g^^'^m^M^i and we find 
B ~ g^^/'^m^M-pi ~ (A/(7)^/^Mpj. Nevertheless, this might not reduce the coefficient B 
significantly since in a models without special symmetry cancellations the constant g can 
not be very large, otherwise loop corrections will induce unacceptably large self-coupling for 
the infiaton. In such models we have g ^ \/\. Moreover, even assuming ^ 1, for the 
effective equivalent temperature which would result to the same strength of the symmetry 
restoration we find T^^ > (A)^/^Mpi ~ IQ-^Mpi. 

The Peccei-Quinn symmetry is restored if i? > (A$/a)/f ~ 10~^^(A<i,/a)Mpi. We see 
that the PQ symmetry is restored at preheating if a > 10^^A$y^, which is a rather weak 
condition. Note also that in chaotic infiationary model the PQ symmetry can be restored 
anyway already during infiationary stage, just due to possible direct coupling of infiaton and 
PQ fields |2^. We conclude that the "thermal" scenario for the axion evolution has strong 
support. 

3) Now let us consider the case of GUT which has large magnitude of the symmetry 
breaking scale $ ~ 10^^ GeV. Since E can not be smaller than m^, the parameter B is 
suppressed for t] particles which have large coupling to $, and which correspondingly have 
large masses Mx- However, at fixed coupling g the decays to instability zones with large 
energy are suppressed, and we possibly can neglect those channels even if the value of g is 
large. In Ref. [jl3[ it was shown that the creation rate of particles does not depend much 
upon g near the surface of zero energy of created particles. Consequently the infiaton will 
preferably decay to particles which have low value of a. Since the particle content of the 
theory is large, we can expect that particles satisfying the condition ~ a^^ < g^^'^m^Mpy 
when there is no additional suppression, can be found. For the case of GUT's this translates 
to the condition a ^ Then, the GUT symmetry is restored if B ^ 10^®(A<i,/a)M|p 
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Using B ~ {\/ gY^'^Mp^ we find the condition a > A$y^. Combining both conditions we see 
that rj particles which will satisfy A$y^ < a < ^ (if exist) can restore the GUT symmetry. 

If parametric resonance is still efficient for particles which does not satisfy the condition 
a < y/g, we obtain B ~ m^M|j/a$^ ~ 10^"^ Mpi/a. Symmetry restoration can be possible 
if A$ < 10-^ 

As final remarks let us discuss the applicability of our approach to the effective potential 
based on Eq. (0). In the usual frameworks developed for equilibrium in Refs. |^ this 



would correspond to one-loop approximation to the effective potential. No significant further 
approximations are made despite the situation is non-equilibrium. As compared to the 
equilibrium, higher loops can be important in the present case, however, and this issue 
deserves separate study. Note also, that the effective potential is not a well defined notion 
in itself, but the effective potential can be defined for the case of special field configurations, 
and this is sufficient for our purposes. The effective potential is equal to minus pressure if one 
restricts to critical bubbles and it is equal to energy density if one considers homogeneous 
field configurations only. Equations of motion for arbitrary field configurations can also be 
easily found in the present approach. It is sufficient to use continuity of the stress-energy 
tensor, dpT^^ = 0, with T^'^ being the sum of Eq. and the stress-energy of the free field 
$. With the use of the Liuville theorem for f{p,x) this continuity condition reduces to the 
equation for the field 

dVo dml f d^p f 
□ $H --{ — - — ^ = 0. (6) 

One of the advantages of our approach to the calculation of the effective potential being 
based on Eq. (or Eq. (^) is that it allows us to find how the coefficient B changes 
when the distribution function evolves according to the kinetic equation |2^, approaching 
an equilibrium. 

In the present discussion we have neglected numerical factors like those which arise due 
to expansion of the Universe. Those will decrease value of B somewhat. On the other hand 
the usually employed description of stimulated decays based on the Mathieu equation takes 



into account only the processes of the form n 2. However, at large phase-space density of 
particles when gf{p) > 1 the processes n m with m > 2 start to dominate. This might 



reduce the evaluation of E of Ref. |T3| and increase the value of B. 

We conclude that physical processes at preheating are very important, especially with 
regard to problems of symmetry restoration, and deserve detailed study. 

When the first version of this work was finished I became aware of Ref. ||2^ where similar 



conclusions were reached. I am grateful to the authors of Ref. for correcting some of 
the statements contained in the original version of my paper. 
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